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Using the method of explicit summation over the intermediate states and 
LS coupling an expression for two-photon absorption cross section in light 
and intermediate atoms in terms of integrals over radial wave functions is 
derived. Two selection rules, one exact and one approximate, are derived. 

In evaluating the radial integrals, for low-lying levels, the Hartree-Fock 
wave functions, and for high-lying levels hydrogenic wave functions obtained 
by the quantum defect method have been used. A relationship between the 
cross section and the oscillator strengths is derived. Cross sections due 
to selected transitions in nitrogen, oxygen, and chlorine are given. The 
expression for the cross section should find usefulness in calculating the 
two-photon absorption in light and intermediate atoms. 


I. INTRODUCTION 


Multiphoton absorption of laser beams by atmospheric atoms and molecules 
and subsequent fluorescence of these particles provide a sensitive tool for 
remote sensing and monitoring of the atmospheric species . 1 To determine the 
strength of the atomic fluorescence, evaluation of the atomic absorption cross 
section is necessary. 

2 

Several methods for evaluation of the cross section are available. A 
suitable method often used is a sum rule or inverse Green's function method 
due to Dalgamo and Lewis 3 which has been applied by several authors ^’ 5,6 to 
two-photon absorption in atoms. The method,, except in resonance absorptions 5 , 
provides reliable results for the cross sections. 

An alternative method is the method of explicit summation over the 
intermediate states. The convergence of the photoabsorption cross section to 
its final value with respect to the addition of intermediate states, as will be 
shown, is rapid. An advantage of the method is that the angular integration of 
the electric dipole matrices appearing in the expression for the cross section 
can be carried out in general for atoms, leading to an expression for the cross 
section in terms of radial integrals. 

Here we use the method of explicit summation over the intermediate 
states. A formula for the cross section in terms of integrals over atomic 
radial wave functions using LS coupling is given. The method is applied to the 
two-photon excitation cross section calculation in nitrogen, oxygen, and chlorine. 

We restrict ourself to cases where the interaction potential between the 
radiation field and absorbing atoms is small compared to the atomic Hamiltonian. 
Then the time dependent perturbation theory can be used to derive the cross 
section. 
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The two commonly used forms of the interacting potential, V(t), are the 

electric dipole form, V(t) =-e£(t).r, with E(t) the electric field vector and 

r the position vector of the electron undergoing transition, and the velocity 

form or the Coulomb gauge form of the potential given by V(t) = (e/mc) A . p, 

where e, ro, c have their usual meanings, A is the vector potential, and p 

7 

is the momentum of the transient electron. 

g 

Lamb has suggested that the electric dipole form should be used, sino„ 

9 

this form agrees with measurements. Bassani et al have shown that for the 
resonance cases, where the difference in energy of atomic levels is a multiple 
of the photon energies and for ls->2s transition, the two forms give the same 
results, but with respect to the inclusion of the intermediate states the 
velocity form converges much slower. 

Kobe has verified Bassani et al results, but has shown that for the 
non-resonance cases the two results are different. However, Grynberg and 
Giacobino' 1 '^ have shown that by taking into account some small terms neglected 
by Kobe the differences between the two forms disappear. 

In addition to the two distinct forms of the cross section due to the 
two forms of the interacting potentials, the cross section has three forms 
depending whether length, velocity, or acceleration forms are used for the 
electric dipole matrices. These forms give identical results if exact atomic 
wave functions are used, but in practice where approximate wave functions are 
used, they can be employed to test the accuracy of the applied approximate wave 
functions. In this way six different forms for the cross section are found. 

For reference these cross sections, except the acceleration forms, are de- 
rived from the Schoedingor equation in Appendix I, treating the radiation 
field classically and the atomic field quantum mechanically. The results are 
identical to results by fully quantum mechanical treatment. Since the 
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acceleration forms are seldom used they will not be given here. Also in 
Appendix 1 an order of magnitude formula is given for atomic N -photon absorption 
cross sections. 

In evaluating the electric dipole matrices a choice should be made about 
the atomic wave functions. It is assumed that LS coupling is valid for light 
and intermediate atoms under consideration. 

However, for transitions involving the ground or the final state and a 
highly excited state the validity of the LS coupling becomes questionable. It 
is more appropriate to use jK coupling^*, where the total angular momentum 
of the core (atom minus the highly excited state electron) is coupled to the 
orbital angular momentum of the highly excited electron, and the resultant 
combines with the spin of this electron. This approximation is not being used 
here. 

For evaluation of the radial integrals, lor low-lying levels single con- 
figuration Hartree-Fock wave functions, and for highly excited states hydro- 
genic wave functions obtained through the quantum defect method have been 
used. It has been verified that the result obtained by this method for the 
excited states agrees within a few percent with the result obtained by the 
Hartree-Fock calculation. Use of the hydrogenic wave functions leads to 
considerable savings in the computer's time. 

The use of the single configuration wave function is not justified for 
a precision calculation. The calculated cross section using this function 
may be off sometimes by factors of two or more from the more accurate multi - 
configuration wave functions. The task of a more accurate calculation is left 
for the future. 

In this paper contribution of the intermediate states that fall in the 

9 12 

continuum have not been taken into account. It has been shown ’ that in 
the electric dipole gauge approximation which will be used here the contribution 
of the continuum intermediate states is small compared to the discrete inter- 
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mediate states. While inclusion of this contribution is straightforward, 
algebraically is cumbersome. It is believed that this contribution should 
be considered only after other approximations which make comparable contributions 
to the cross section have been removed. 

In Section II that follows angular integration over the electric dipole 
matrices is carried out, and the cross section is expressed in terms of the 
radial integrals. Two selection rule.' ...re also derived, and a crude method 
for cross section calculation using ';.he oscillator strength is given. In 
Section III cross sections for selected transitions in N, 0, and Cl are given. 

For the specific problem of the two-photon excitation in atomic oxygen, in 
Appendix II population of the excited and ionized states of the excited atoms as 
functions of the laser pulse duration, the collision time of the atom with the 
ambient gas in the atomosphere, and the atmospheric opacity to the probling laser 
beam have been calculated. Results are presented graphically and through a 
table. 

II. FORMULATION 

A. Angular Integration of the Amplitude 
We carry out the angular integration over the amplitudes given in the 
Appendix by Eq. (A14) . Since we are dealing with light and intermediate atoms, 
we assume that LS coupling provides a valid coupling scheme. The state la> is 
then given by 

|a> = | Ca 1 S 1 L 1 n£)SLJM J > Cl) 

In this expression n£ is the principal and angular momentum quantum numbers 
of the active electron undergoing transition, is the spin and orbital 

angular momentum quantum numbers of the residual atom, and otj represents all 
other quantum numbers of the residual atom. SLJMj are the spin, orbital, and 
total angular momenta of the atom, and Mj is the projection quantum number of J. 

The evaluation of the dipole matrices can easily be accomplished by the 
use of the powerful method of Racah. In evaluating the matrices in (A14) we 
make use of (1) to expreso |i>, |m>, and |f>. We also use the notation 
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E.e.r = P (10) where r. are the position vectors of the atomic electrons, and 
1 i 

P (10) is the z-component of P = L r (Ref. 13) . It should be noted that 

i 

summation over m stands for a summation over the principal, azimuthal and 
magnetic quantum numbers of the intermediate states as well as an average with 
respect to the initial, and a sum with respect to the final magnetic quantum 
numbers . 

We designate the initial, intermediate, and final states by the unprimed, 
double primed, and primed symbols, respectively. Then (A14) can be written 


j* _ / 


♦ * , / ,-a- 


(*) j x I J — iSLJSL Jsl jj / ^ 

= fn yE - £ / 

where T is defined in the following. 

. 14 

For the case of non-equivalent electrons in a shell T is given by 

faWK'j'; L, J5LTS LJ S LJ') 


= -L- Z (l I He. IHe. r/.> 

= -L- 7Z M <(4,s,L, -»$)SL iMylroojj^s, Lr"/j s lj\„} 

* <(K s , L r*fJ s l" J \ l ^ (l0) l ( *‘ s l L r' J> J SLJ 

^s-tL 4 L - +i '(z j+ij ] A 

[£ (v, '•>)(» ■» 

A7 7 3 7 3 ' J 
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( 3 ) 


All symbols have their usual meanings and 

r°° 3 

<*,/hpo)ii«v > - J \ l ; j V;/ r 

o y 


where Z > is the larger of Z and Z , and and R^, are single electron 
radial wave functions. 

In the above we have assumed that the residual atom designated by the 
quantum numbers is "frozen", and does not change its state during thr 

transition. Similarly, the total spin angular momentum S does not change 
during the transition. 

If there were n equivalent electrons in a shell, the expression for T 
14 

becomes 



/ 


( 5 ) 


where T is given by (3) and the arguments of T and T are the 

same. The bracketed term is the coefficient of fractional parentage in which 

£ n a SL are the quantum numbers of the shell that has equivalent electrons, 
n ^ 

and Z (o^S^Lp ZSL are the quantum numbers of a hypothetical shell with n 

electrons in which there is an active electron with angular momentum quantum 

1 

number Z, a residual atom with quantum numbers (ot-jS^Lj) , and the shell 
has the total spin and orbital angular momenta quantum numbers S and L. 

Tables of numerical values of the coefficients of fractional parentage are 
given by Racah^ and Rohrlich"'^. 

We can combine (3) and (4) with (21 for the case of non- equivalent 
electrons, and (31 , (4) , and (5) with (2) for the case of equivalent 
electrons to express the cross section in terms of rad-' al integrals. 

I ! It II II 

In addition, sincr variation of E(n Z L J ) with respect to J' * for 
light and intermediate atoms is slight and of the order of the fine 


II * * If II 

structure constant, we take the average of E(n £ L J ) with respect to 
J'', and sum only the numerator in (2) with respect to J M . Then the cross 
section can be expressed as 


CT(SLJ-* 5 L'j) -?n YE 


! I d i/cfR h^/jR(~yy/jj 
-TV / 


a. 




( 6 ) 



(7) 


where C(£ M L") are constants expressed in terms of 3-j, 6-j symbols, and 
the coefficients of fractional parentage, and are obtained by comparing f6) 
and (7) with (2) through (5). For the case of non-equivalent, electrons it 
is given by 


d(/L*) = H 


JJ+/+ 2.S+L+L- 




y i 





7 7* I 

«l 7 -* T 





(8) 


i i iii 

where £^ is the larger of 1 and £ , and Z^ is the larger of £ and £ . 

For the case of equivalent electrons anologous to Eq. (5) the right-hand 
side of (8) should be multiplied by the coefficient of the fractional 
parentage. 
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B. Selection Rules 


1. An exact selectio n rule . We consider the express on in squared 


bracket on the right-hand side of (3) consisting of a summation with respect 

to Mj. Since Mj is a demi -variable, the expression does not change in value 

1 7 

if we let Making use of the properties of the 3-j symbols it 

follows that 

j-H j" i j / - r 7 " I I j-fy a" / j/jVf ) 

M \/M -M oJ\M -ft o) M l-AI o]\-m /Mj o) C c 

J J J j J 


+ j— f ' 3 J " 1 j/ 3 ' J / 
6-[ *L-M <,)(#; -« r o 


Therefore the expression vanishes unless J+J'+2J' ' is an even integer. For 
J, J', and J' 1 integers this implies that J'-J should be even. Then of all 
permissible values of J' only J' = J, J+2 are allowed. 

For J, J', and J' ' half odd integers, the condition J+J'+2J'' even 
also implies that J'-J must be even. Then again J'=J, J^2. We therefore 
have the following selection rule: 

J ! - J = 0. + 2 (1C 

This rule is in addition to the usual selection rules for single transitions 
in which A£ = +1, AL = 0, +1, and A.T = 0, +■ J , where l, L, and J have their 
usual meanings. 

2. A n approximate sel e ction rule . This approximate selection rule 
arises when we neglect splitting of a level with a given S and L due to 
different J values, i.e. , the fine structure splitting, and applies to the 
case when J’ = J. With the neglect tf the J splitting the denominator inside 
the squared bracket on the right hand :• ide of (2) becomes independent of J'', 


and the sum with respect to J' ' in (.2) car. he cone analytically. 


For J' * J the sum with respect to Mj on the right-hand side of (3) 
becomes equal to 1/3 (l^f. 17). Then the sum with respect to J*' in (2), 
using (3) , becomes 




ztt+l) 

— (■“ ) (iL^l 



17 

Eq. (11) was derived using the ortho gonaliz at ion of the 6-j symbols. 

Eq. (11) implies that 

L» -L = 0 when J* - J = 0 (12) 

i.e., when the initial and final total angular momenta are the same, the 
initial and final orbital angular momenta should also be the same. It will be 
shown in the next section that transitions for which J'=J have large cross 
sections. When (12) holds, (3) reduces to 


i ■" / // / ♦ it t / . 

iC-v. L, iSLJSLTSLJ/ 


Using (2) we see that cross sections due to transitions which are 

4,22 

forbidden according to (12) are of the order (Z /y ) times smaller than 
allowed transitions, where Z is the effective charge acting on the transient 
electron, and y is the fine structure constant. 


(13) 


C. Method of the Oscillator Strength to 
Calculate the Cross Section 

18 

A method used by Mcllrath, Hudson, Aikin and Wilkerson employs the 
available values of the oscillator strengths to calculate multiphoton 


- 10 - 




excitation in atoms and molecules. The approximation in the method consists 
in neglecting all intermediate terms represented by the summation with respect 
to m in (A14) except one term with the largest oscillator strength. In this way 
an order of magnitude cross section is calculated. The advantage of the method 
is that no elaborate calculation is necessary, and use is made of the most 
accurate oscillator strength values available in the literature. However, 
when no strong line occurs in the transitions considered, the calculated values 
may be fur off from the true values. Examples will be given for N and 0 using 
this method of calculation. 

19 

The oscillator strength between two levels Oj and ctj, is defined by 


X j— 2ZI<(<S,L r /)SL7M J IZJ. rl«S,Lj/jsL'7rtSf 

a J* 1 Mm * ' 

where AK(cCj, u^,) is the energy difference between the two levels, and 
other terms have been defined previously. Eq. (14) can be evaluated from (3) 
by letting in this equation the primed quantities to become unprimed, and the 
double primed quantities to become primed, and neglecting the phase factor: 

- y &E(<* Jj °(j,)(2Li-l)(2.L / -H)( J - 

When equivalent electrons are in a shell, the oscillator strength is 
obtained by multiplying the right hand side of (15) by the factor 


n(£ n aSL 


z' 1 (o:^S ,I.p iSL) ‘ , all notations being defined following eq. (5) . 


We now introduce tin line strength 9c -i } and the multiplit strength 

16 


tfCM) by tlx* following equations 10 : 


(14) 


¥( 


(*T>n(2.j+tj f 5 j 




*1 


5 t"/ 


I 

L~ 



(16) 




(17) 


¥(cM)~as+o(iL+o&L'+oJi(*£ofrL 

When equivalent electrons are in a shell, the multiplit strength is obtained 

by multiplying the right hand side of (17) by the factor n(£ n otSL j£ Jt n * (o^S^L^JiSL) 2 . 

13 

The absolute line strength is now defined by 

= $?U ) <- fn / (18) 

The line strength is normalized such that 

21 tf(jL) = / (i9) 

JJ' 

It the;; .follows that the absolute multiplit strength is given by 

s c7+j% (20) 

Jj' 

Through Eqs. (15) -(18) the oscillator strength for transitions between two 

20 

levels could be expressed in terms of the absolute line strength 

/fay** t21) 

An approximate expression could also be obtained for the oscillator strengths 
tor transitions between two multiplets. The approximation consists in 
neglecting the energy differences of different levels within a multiplet. 

Using the definition 

/ K — (‘iS’+OfaL+O 2Z( 1 J+0 f (22) 

v LJ JJ, 

and Eqs. (20) -(22) we find that 

AE ^L^t:) t23) 
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where An (a L ,a L ') is the average value of AE(Oj,aj,) with respect to J and J*. 

If only one intermediate term is included in the sum in (A14) , this 
equation can be written - 


j'L . 

* ttj ' 


J 

2 -T+/ 


/ZT^f 'lP(i»Jir*xn?«ojir>l 

A 1 M 


(24) 


where we have averaged with respect to the initial, and summed with respect t - " 

t I ! 

the intermediate and final magnetic quantum numbers, and T, F and F stand 
for (cijSjLjnajSLJHj, (o^S^n’ Y ')Sl' ' j"Mj" and (ct^I^n Y)Sl/ j'm j( , 
respectively. 

Eq. (2/n can be expressed in terms of the oscillator strengths by 
making use of (3) and (15): 


(i) <ivY je z 


if 

2-1 tf 
2-J+l 


{j*±2^d 

A£ A£ ^j-^ < rv' J 


7 T* / \/ 7 

A1 -/»L 0 j V M ~ 

J J ' 3 

In the above £ is a number of order unity. For .J=J', r, =1. Eq. (25) is the 
desired equation for calculation of the cross section. 

21 22 

Values of railtiplet strengths have been tabulated by Goldberg 5 . It 
is interesting to note that the tables given by Goldberg are contained in the 
simple expression due basically to Racah given by (17). In the tables, however, 
values of the mult ip let strengths are not listed for separate parents for the 
case of equivalent electrons in a shell, rather values corresponding to sum 
over the parents are given, while i-q. (17) giver, vaiue of the multi plot 




strength for individual parent. Also, tables by Goldberg give the square of 


the transition amplitudes, while in the present calculation values of the 
transition amplitudes are necessary. The use of the Racah's formulas for 
multiphoton absorption cross section calculation then is a necessity. 

III. RESULTS AND DISCUSSIONS 

Formulas derived in the previous sections have been applied to the two- 
photon excitations in atomic nitrogen, oxygen, and chlorine. In using the 
explicit summation method, for low lying intermediate states use has been made 
of the single configuration numerical Hartree-Fock wave functions^' . For 
high-lying levels we use hydrogenic wave functions with non-integer principal 
quantum numbers fixed by the quantum defect method. This is justified since 
for highly excited atomic states hydrogenic description of the states are 
valid. 

The hydrogenic wave functions with non-integer principal quantum numbers 

used here is in the form of an asymptotic expansion discovered by Eddington 

20 

and Suguira as reported by Bates and Damgaard. This function, not being a 
hydrogenic eigenfunction, diverges at the origin. However, for exact radial 
wave functions the value of the integrand in the integral given by (4) vanishes 
at the origin, and for highly excited states has small values for distances 
from the origin of the order of the Bohr radius a Q . Then to circumvent the 
problem of divergence, we introduce a lower limit cutoff in the integral (4) 
equal to a Q . It can be shown numerically that the value of the integral is 
insensitive to the slight variation of this cutoff. 

In evaluating the radial integrals, states are divided into two groups, 
according to whether the values of their wave functions are given through 
Hartree-Fock, or hydrogenic wave functions. It is tested that values of 
the integral does not change appreciably by switching some states at the 
boundary of the two groups from one group to the other. 


We choose transitions which lead io strong emission lines, since fluore- 
scence from these lines are being detected. The initial state as a rule is the 
ground state. The final state is chosen depending on the range of the laser 


frequencies and allowance of the selection rules. In this way only a few 
suitable transitions can be found for each atom. 

3 4 o 

For the case of atomic nitrogen we consider the two transitions 2p S 
-*■ 2p 2 ( 3 P)3p 4 *°yo, and 2p 3 4 S° 3/ , ■* 2p 2 C‘ i P)3p 4 n° 7/2 . Values of J' in the 
second transition could be 1/2, 3/2, S/2, 7/2. However, J' = 1/2 and 5/2 are 
forbidden by the selection rule (10). Also J' = 3/2 is forbidden by the 
approximate selection rule (12). Then we are left with J' = 7/2 only. The 

t t IV 

coefficient C(£ i ) from ( 8 ) for the first transition are C(01) =-0.19245 and 
C(21) =-0.35490, and for the second transition are C(01) = 0.14903 and C(20) = 
0.029805. 

Table 1 gives values of cross sections for the two transitions. For each 
case 8 intermediate states have been employed, and contribution of each inter- 
mediate state to the total cross section is shown in the table. Addition of 
more intermediate p* aces would change the cross section in the third significant 

4 

figure. Hxpressee m cm , we find that for S,, , 0 transition the cross 

-35 4 

sect ion is 2.74x10 ‘ cm . while for »S„ ... - D - /0 transition the cross section 

>V 2 //2 

is 3.97x10 ' cin , almost a factor of 7 smaller than the first transition. In 

ill cases considered it seems that the cross section is largest when the 
initial and final orbital and total angular momenta are the same. 

Using the method of the oscillator strength, cross sections for the two 

• 30 — 3(3 *4 

transitions by Fq. (23 l are 5. 2x1 J ‘ and 6.9x10 cm . 

4 3 4 o 3 

for the case of at onuc oxygen the transit lens 2 p **P,, 2p ( ) 3p * P 

- o , - 

have neon considered. Of the throe values, of J ! ~0, 1, 2, the value J'=l is 
excluded through chc select ten rule (10). Using (8) the coefficients C(V * L* * ) 




T 


‘ 


for J'=0 are C(00) =-0.036289 and C(22)=0. 1283, and for J'=2 are C(00) =-0.0072578 
and C (22) =0.2566. 

Table II gives cross sectional values in oxygen for the final angular 
momenta J'=0 and 2. Similar to Table I, eight intermediate states and their 

4 

contributions are included. Expressed in cm , the cross sections 0 q (O), 

a Q (2) and a Q ^, are given by 1.71x10 ^ , 6.60x10 and 6.77x10 ^cm^, 

where a „ is the sum of the two cross sections. It is seen that the cross 
ol 

section for the transition in which the total angular momentum does not change 

is 39 times larger than when it changes by 2 units. 

£ 

Pindzola has calculated cross sections for the above transitions in 
atomic oxygen, using the inverse Green's function method of Dalgarno and 

3 

Lewis. Using the electric dipole form of the interacting potential he obtains 
-36 4 

a value of 7.62x10 cm for Considering the many approximations used in 

the present paper, the agreement between the two independent methods, if not 
accidental, is striking, and is an indication of the reliability of the cross 
section. 

These cross sections can also be calculated using the method of oscillator 

24 

strength. Using the value of the oscillator strength for the single inter- 
mediate state 3p^( 4 S°)3s^S° and Eq. (25) we find for CT^ (0) , (2) , and 0 qT the 

_ <7 7 _ •7/1 _ rr/L A 

values 1.88x10 , 2.35x10 , and 2.54x10 an respectively. Therefore the 

value of 0 qT by this method is by a factor of 2.7 smaller than the value of 
0 o> j, obtained by the method of explicit summation. 

The value of 0 q ^, using the method of oscillator strength given by 
Mcllrath et al^ for 0^ is 7. 13xl0 _ ^^cm^ which is by a factor of 2.81 larger 
than 0 q ^, given here. The discrepancy is due to the neglect of some weighting 
factors in (25) by these authors. 
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i • _52 _o 

In the case of atomic chlorine the transition considered is 3p P.^ 
3p 4 ( 3 P)4p 2 P° /2 Using (8) the coefficients C(&' 'L' ') for the 3 intermediate 
series are C(011=-l/9, C(21) =-1/18, and C(22)=-l/6. Table III lists 11 


intermediate states and their contribution to the total cross section. In 

4 -36 4 

terms of cm the total cross section is 1.89x10 cm . 


IV. SUMMARY AND CONCLUSIONS 

Using Racah's formulas an expression for two-photon absorption cross 
section for light and intermediate atoms in terms of integrals over radial wave 
functions is derived. An exact and an approximate selection rules, in addition 
to known selection rules, are found. The formulation is applied to selected 
transitions in atomic nitrogen, oxygen, and chlorine. 

When oscillator strengths calculated using single configuration Hartree- 
Fock approximation are tested against more accurate tabulated values, on 
occasions they are different by factors as large as 2-3. This suggests that 
cross sections calculated here may also be different from their accurate values 
by these factors. 

In a more accurate calculation in the future multiconfiguration Hartree- 
Fock wave functions should be employed, jK coupling scheme should be considered 
for the excited states, and continuum intermediate states should be included. 

In addition, length and velocity forms of the cross sections as given in 
Appendix I should be employed to test the accuracy of the wave functions 
employed. 
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APPENDIX I: DEVIATION OF TWO PHOTON ABSORPTION CROSS SECTION 


The transition probability for a particle to go from an initial state 

li> to a final state lf> when the interaction potential is turned on from time 

27 

t to time t is given by 

W(i-f) = |<f |u(t,t Q ) | L>| 2 (Al) 


where U(t, t Q ) is the evolution operator. The transition amplitude in 

different orders of the interaction potential can be expanded in the following 
27 


way 


oo 

<fiuct, t n )|j> = y~ <fiu (n) (t, t o )ii! 

n=l 


(A2) 


The two-photon absorption process corresponds to the second order amplitude 
.27 

■* t, 


given by 


(// U% 7Z [/*, fe'h IV“* 


(A3) 


where we have introduced = E , E^ being the eigenvalue of the eigenstate 
lk>, and V„ (t) are the expectation values of the time dependent interaction 
potential V(t). 

In empty space the vector potential of the radiation field is given by 


a ~ a 0 (* • r [-* a • r — w tjj 
= A» [e~+ ■* e ™-1> (<'“*;] (M 

where k is the propagation vector, and w is the frequency of the field. The 
approximate equality in (A4) results from the fact that in the visible range 


of the spectrum k»r«l. 
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The first terra on the right hand side of (A4) corresponds to the absorption 


of photons from the field, and the second to the emission of photons into the 
field. 

Since we are dealing with absorption only, the first terra will be used in 
the calculation. The interaction potential in the length form is then given by 


Vl-t)-=-e 7 ZE(i)' r. 

*• 4 . <3 


V 0 - \ ZT/4 0 ' r d J (A5) 

where is the position vector of the ith atomic electron. Substituting 
the potential given by (AS) in (A3) and making the plausible assumption that 
t -»«> we find that 

<flU w h> =- * )> 



~°Q -oO 


'Trj "k C H’ ~ ^ ~+ c £ ~ ~t X CJ -f 2 t gj 

* * f 


/ 


(A6) 


where we have introduced the converging factor e. At the end we let e-K). 

The choice t •*•••” while convenient is not valid when the perturbative time is 

“17 2 

comparable to the periods of atomic, orbitals (Of the order 10 / Z" sec, where 

Z is the effective charge acting on the transient, electron) . 

Let us assume that the first order transition is zero. Then by letting 
2e->-e we obtain from (Al) and (A 6) 





/r 


x 

/ Pl\/ \ \ / /./ / i *€? 

\ q » i'^y \'*» i V<>l i\ / £. 

^ ~ - * £./ i j ( Cd , - L> -f zu]\. g x 

1 f 
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Using the definition of the transition rate and by letting e-*o we obtain 





(O . tO 60 + t £ 



(A8) 



We now substitute the value of V Q from (A5) , but first we relate A q to the 
average field intensity T which is an observable quantity, while A q is not. 

From the definition of the Poynting vector we have that 



EkH 


4. 

CO 

2. 7TC 
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Making use of (A5) and (A10) in (A9) leads to 

A 


K - . H P 1 to.- -f. to +<■'€ ' * 

% Orj i T 


e being a unit vector in the direction of the polarization vector A q . 

The transition rate is related to the two-photon absorption generalized 

(2) . 

cross section a by 


)?= <r ,v (i/t») 


Z 

(A12) 


(All) 


Combination of (All) and (A12) leads to 


2 . 


(A14) 







/ 


where y is the fine structure constant and E =t w is the energy of the photons. 

f21 ( 2 ) 

We notice that a Q K J has the dimension of 4th power of length while CJ 

has the dimension of 4th power of length multiplied by time, and R has the 

dimension of inverse of time. 

The converging factor e introduced in (A6) has a physical meaning. It is 
proportional to the width of the final state. The 5-function in (A13) can then 
be shown to be approximately equal to 2x f /iT, where is the lifetime of the 
final state. ^ 

Equations (A13) and (A14) are identical to similar equations derived by 

26 

fully quantum mechanical methods. 

Making use of the commutator relationship 

e - ~ (Z o (£ I e ’ r I (Ai5 

where E n = E - E, is the energy difference between states la> and lb> 
ab a b ° 

in rydberg units, and a^ is the Bohr's radius, (A14) can be written in the 
alternative form 



Energies in (A16) are expressed :u rydberg. While for exact atomic wave 
functions the two expressions i and (A161 are equivalent, they differ for 
approximate wave functions. The degree of the difference is usually taken as 
a measure of the accuracy of the approximate wave functions used. 


mi .. 


If we had to use the velocity form or Coulomb gauge form of the potential 


given by 


V(t) = (e/mc) 


(A17) 


we would obtain the following value for a 


( 2 ). 


3 i. 

(T =?n(a 


h , /^_ qr^_ y.hx-iK*- */«> i 

r /' /£ jZ- E.-E^+E r 


2. 

(A18) 


Making use of the commutator relationship (A; _J , (A18) can be written 



(A19) 


We therefore have four equivalent forms for the cross section. These 
forms can be used to test the accuracy of the wave functions used. 

According to the principle of microreversibility if the total Hamiltonian 
remain invariant under the time reversal, the transition probability in the 
reverse direction is the same as the direct transition. 

H(-t) = H(t) R(f,i) = R(i,f) (A20) 

It is seen from (A4) that the time dependent part of the Hamiltonian is 
invariant under the time reversal. The reverse transition corresponds to 
the emission into the field, and the creation operator expfiwt) should be 
used. The cross section for the reverse process is obtained by letting in (A14) 
i-t-f and ur>--a). Then (A14) is transformed into 


<T 


r } = fv/V/T: ^ ^ / Wi) 

p *>1 
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Eq. (A21) is however not an independent equation, since by the conservation of 


energy E--E -E => E.-E +E, and (A21) is the same as (A14) . 

± m 1 m 

Making use of the closure relationship, an upper bound and an order of 

(2) 

magnitude value can be found for a Q . According to this relationship 

2IY/ le.? I™) (■»> l € • r*/ *■> = Y I t^l ‘ > (A22A) 

where |i> represents the total wave function of all electrons in the atom or 
the molecule. However, when, one electron wave functions are used, as usually 
is the case, (22A) should be replaced by 

ZT <//«. 

'hf V ' J (A22B) 

- 2 Kf IS. ?!*><» i* 

M 

where M stands for occopied atomic or molecular shells inaccessible to the 
transient electron. Eq. (A14) can then be written 


(£ (i! 7(iffYa?} [ 


(E, - a ^ £ L- 
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'tot-Y, 
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(A23) 


where in (E^ -E -E) m ^ n , E m corresponds to an intermediate eigenvalue which is 
closest to 1/2 (E^-Ep . 

Since the energy ratio and terms inside the bars in (A23) are of the 
order unity, the two-photon absorption cross sections in atoms and molecules 

2 n n, _35 4 

are expected to be of the order of 2ir(27TYa o 1.035x10 cm . For the 
three atoms N, 0, and Cl considered hero, indeed the calculated cross sections 
differ from the above number within a factor of 5. The energy denominator 


in (A14) as a rule does not vanish, since E m for any m can only accidentally 
fall halfway between E^ and E^. 

The above order of magniti’.de calculation can be generalized to N-photon 

2 

absorption cross section calculation. Using the formula for this cross section 
the order of magnitude expression will be given by 



in (itiY cl* ) "(t At A 

in (l.zfmo 1 *) (h. In *!o l7 j 
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Thus for the three and four photon absorption cross sections we get respectively 
6.41xl0~^°cm 8 sec and 3.97x10" ^cm 8 sec 2 (Compare Table I of Ref. 2). On the 

other hand the single photoabsorption cross section is of the order of 

. 2 2 0 n . . _-18 2 
4 tt Y a = 8.04x10 cm . 

1 o 

Since the derivation in this section was based on perturbation theory, 

it is in order to find the validity of the theory in dealing with intensities 

of short pulsed laser beams. The validity criterion is that eE Q <r> a «V a , where 

E q is amplitude of the electric field of the laser beam, <r> a is the average 

atomic radius, and V is the atomic potential. From (A5) and (A10) we find 

that |E 0 | =)M7 c. We also assume that <r> a is of the order of the Bohr 

radius a Q . If we then express the energy of the laser beam e in mJ, its 

2 

cross sectional area A in cm , and its pulse duration in n sec, the value 
of [ E q j <r> a in volt will be given by 


l£. I < r > a (y) = *• 7 ® VjtS 


€ (<*7 Tj 


(A25) 




Since the atomic potential is of the order of electron-volt, for the perturbation 

theory to be valid we must have |E Q |<r> a «l. As an example, for a typical 

-4 2 

case of e = 1 mJ, A = 10 cm , and = 10 n sec, we find the laser interacting 
potential to be 0.00229 eV, indicating the validity of the perturbation theory. 
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APPENDIX II: POPULATION DYNAMICS 


A. Rato Equations . With the knowledge of cross sections for different 
processes, we can calculate population of different states at the end of 
a laser pulse. For simplicity we assume a three level atom consisting of 
the ground, excited, and ionized states. The processes that should be taken 
into account in the rate equations are the following: (a) Resonance two- 

photon absorption by the ground state, (b) Resonance two-photon emission by 
the excited state, (c) Spontaneous emission of the excited state, (d) Photo- 
ionization of the excited state, (e) Triple and higher multiphoton coherent 
ionization of the ground state, (f) Radiative recombination of electrons and 
ions which follows ionizations, and finally (gl Population and depopulation 
of atomic states through collisions with the ambient gas. 

It will be shown in Section B that for laser duration of interest (e) 
can be neglected. Similarly, radiative recombination lifetimes are too long 
compared to the short-pulsed laser beams, and (f) can be neglected. In 
dealing with more complex atoms than those dealt with here a fast recombination 
process, namely dielectronic recombination, occurs which should be taken into 
account. In this case autoionizing processes should also be accounted for. 

Here we neglect these processes. 

Finally it will be shown in Section C that for gas densities less than 
10 1 /cm'" 1 corresponding to the atmospheric conditions above 30 Km the collision 
effects can be neglected. 

Additional assumptions are that the two-photon excitation is the only 
mode of excitation of the ground state, and ail electrons that cascade by 
the spontaneous emissions end up in the ground state, i.e., no metastablo 


states exist. 


Kf 


We are then left with four processes (a) through (d) . Let Nj , N 2 , N 3 
represent the instantaneous population of the ground, excited, and ionized 



states, and N q the initial population of the ground state. Then Nj+N 2 +N 3 * N q . 
The rate equations for these states are given by 

0)hl, /(ft — - ■+ } ~ ~ ^ir 

d</di - < ( v, - 4^ a/ l , (fm) CA26) 

W/di = ^ - R lr -r + (f -) , K 

r o’) 2 

In these equations R 2r = TgI N is the rate for two-photon absorption 

or emission. Tg is the smaller of the inverse of the laser's band width or 

18 

the atomic Doppler width , and I N is the number intensity of the laser's 

n _ 1 

beam in units of cnT^sec . Rj = cr^ I N is the rate of single photoionization 

of the upper state with cr^ the photoionization cross section, and finally 

A(f,m) are the Einstein's transition probabilities for going from the upper 

state f to the intermediate states m, where E <E~. 

m x 

The initial conditions at t=0 are that Nj=N q , N 2 =N 3 =0. Strictly speaking 
the values of N 2 and Nj at t=0 are determined by the Boltzman factors. How- 
ever, for conditions of interest in the laboratory and stratosphere the error 
due to the neglect of these factors is extremely small. 


Subject to the above initial conditions (A26) can be solved by the 
standard methods leading to 



x being the duration of the laser's pulse. From the above it is easily 

verified that N.+N 0 +N_ = N . 

i l 3 o 

Since the validity of (A27) is based on the validity of the perturbation 
theory, x or the duration of the laser's beam must be much larger than the 
orbital period of the running atomic electron (Cf. Appendix I). 

The above equations are applied to find the state population of atomic 
oxygen due to interaction with a laser beam of energy 1 Joule, a radius of 
the cross sectional area at the position of interaction of 18 cm, and a 
variable pulse length. 

( 2 ~\ 

The parameters needed to evaluate (A27) are cr Q , cr, and EA(f ,m) . Values 
f 21 

of a J are given in Table II. Since is the photoionization cross section 

25 

of the excited state, the quantum defect method employing hydrogenic wave 

functions can be used to find its values. 

In calculating we notice that initially the electron is in the state 

2p^( 4 S 0 )3p^P. The final states then should be 2p 8 ( 4 S°) es^S° and 2p^(^S°)ed^D°, 

24 

where e is the energy of ejected electrons. Using tables of energy levels 

we find by extrapolation that the quantum defects of es and ed arc approximately 

1 and 0. Then the photoionization cross sections for these two states and 

-19 -19 

their sum can be calculated and are given by 0.972x10 , 0.119x10 , and 

0.109xl0~ 18 cm 2 . 

Using the values of cr ^ and cr^ we find that = 8.45x10 ^Xg/x 2 , and 
Rj = 3.16 x 10~ 4 / t . 

In calculating E m A(f,m), with good approximation we can replace this sum 

3 3 3 3 

by the transition probability for the strong line 2p 3p P 2p‘ 3s S with an 
oscillator strength of 0.898 (Ref* 24), leading to a transition probability 
of 8 . 42xl0 2 /sec. 

With R„ , R- and E Aff,m) so determined, (A27) can be evaluated for the 
2r 1 m 

cases x t1 =x, and x =x„, x n being the inverse of the Doppler width. Results 

o o u IJ 

are presented in Figures 1 and 2. 


In Figure 1 the form R 2 =o ^ 2 \l N 2 bas been used for tbe two- P hoton rate ' 
As is seen both N 2 /N q and N,/N n have constant values for very short pulses up 


3' o 


,-10 


to t= 10 sec, then they fall off exponentially. In the regions that these 

-7 

ratios are constants, (A27) take a simple form. For t« 10 , it can be 

shown that A ± t« 1, and (A27) reduces to 

fA a I - KrT[i- i (f+A- K,) r] 

(A28) 


R ir T [l-i (A^A'jT 7 

Vj/zv; ^ \ 


showing that N 2 /N q and Ng/N^ approach constant values for short pulses. 
The expression for N 2 /N q given by (A28) ± s in full agreement with a 


18 


similar expression given by Mcllrath et al , since their expression is the 
limiting value of N 2 /N q given by (A27) . 

( 2 ) 2 

In Figure 2 the alternative form R 2r =a 0 t d N n has been used * In this 

-15 

figure three distinct regions are noted. For T<10 both ratios are constants. 
For 10 -14 <t< 10 -9 the ratios fall as the inverse of x, and for x>10 7 the ratios 

9 

decrease as 1/x . These behaviors can be seen from (A27) as follows: 

When 'f g =T D , R 2r =1.07xl0 -15 /x 2 , while R 1 and E m A(f,m) have the same 
values as before. The following cases follow: 

(i) x<10~^. Eq. (A27) shows that A + x»l and A x«l, leading to the 
following limiting forms: 




(A29) 


Then to first order the three ratios have constant values in agreement with 
Figure 2. 
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(ii) 10 <x<10 . In this case A T«l, leading to the following 

limiting forms: 

N, /N 0 m I y K/K « R ir T , Sr i- R ir K , T (A30) 

Then N 2 /N q and N 3 /N q fall as 1/x in agreement with Figure 2. 

_7 

(iii) t> 10 In this case £A(f ,m) >>R^>>R 2r , leading to the following 
limiting forms: 

Af/V ci./ a£ /v, R ir /TIM fa) 

3 /__ wu 

2 

Then N 2 /N q and N 3 /N q fall as 1/x in agreement with Figure 2. 

18 

Case (ii) for N 2 /N q is in agreement with the result of Mcllrath et al. 

-14 -9 

Therefore their result is valid for 10 <x<10 , which is the important 

range of x in application. 

The results for N„/N and N_/N given in Figures 1 and 2 are combined 

tL 0 o 0 

in Figures 3 and 4. It is seen that the ratios N 2 /N q and N 3 /N q have almost 
constant values for x<l . 266x10 sec. Then they fall off rapidly. The 
constant values are 0.845x10"^ for N 2 /N and 1.335x10 ^ for N~/N o . 

B . Effect of Three - Photon Ionizations 

The energy of the laser beam is sufficient to ionize the ground state 
through three-photon ionization. This process is different from the two- 
photon excitation of the ground state and the photoion.ization of the excited 
state in two respects. The former is a coherent, while the latter is a 
non-coherent process . In addition the former process depends on the 
population of the ground state, while the latter depends on population of 
both the ground and the excited states . 
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When the three-photon ionization of the ground state is included, the 


last equation in (A26) should read 

d V dt = r x n 2 * r 3 n x 


(A32J 


where is the rate of three-photon ionization. We substitute and ^ 
from (A27) into (A32) and carry out an integration with respect to t. The 
following result follows: 


r-u o-* A ~v] 

K A ' A s T) * A t / s Tl )W 

+*,[ !S |J7 S 

For short pulses, similar to the derivation of (A28) , A T<<1, and (A33) 
reduces to 




(A34) 


C3-) 3 

Realizing that Rj=a^ •'xl^ , the ratio of the second term to the first on 

the right-hand side of (A34) is given by 2o / a a . Making use of (A24) 

2^ -18 2 

and the fact that a ^ is of the order of 4irya o ^2.56x10 cm , the value of 

-17 

this ratio is given by 4.84x10 /x. 

-17 

We therefore conclude that for x»10 sec, the three-photon ionization 
can be neglected. Similarly, higher order processes can be neglected. 


C. Collision Time . Collisions of molecules under study with the ambient 

gas have dual effects. They broaden the absorption and emission lines, 

and have quenching effects on the excited states . 

To obtain a rough estimate of the collision time we use the classical 

picture that molecules are made up of hard spheres. Then the collision rate 
28 

is given by 



_____ i ___ 2. t'/* 2 - 

/T 77-^ (r, + n_J (-ir -+ VI J 


(A35) 
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where r, and v. are the radius and average velocity of the molecule under 
study, and n 2 , r 2 , and v 2 are the density, radius and the average velocity 
of the ambient gas. 

We use (A35) to calculate the collision rate of the excited state 

2p^(^S°)3p^P of oxygen with the ambient gas. The average value of the radius 

of the excited state of oxygen according to the Hartree-Fock calculation is 
-8 

given by 3.67x10 cm. For simplicity we assume that the ambient gas is ma^e 

up of 80 percent of N 2 and 20 percent of 0 2> Based on the hard sphere model 

and the viscosity data, and neglecting a slow varying temperature factor, 

—8 “8 

the radii of N 2 and 0 2 are given respectively by 1.88x10 cm and 1.81x10 cm 
(Ref. 28). Using (A35), the rate coefficient k, and the collision time t £ 
are given by 

0.^37 x'o p-j-’ 

) (A36) 

I ^ -=//0 ~ 2 . 2-7 mo 

29 

Using the atmospheric data for densities and temperatures , (A36) have 
been evaluated. Results for as functions of altitude are given in Table IV 
and Figure 5. It is seen that in the range of 20-120 Km the collision time 
is larger than the Doppler time 'v 1.266x10 ^sec. However, for altitudes 
less than 30 Km, the collision time becomes comparable to the average lifetime 
of the excited states. For these altitudes the rate equations in Section A 
should be modified. 

D. Opacity Calculation 

It is of interest to estimate the absorption of the probing laser beam 
before it reaches its target. The main absorbing processes are the resonant 
two photon absorption by atomic oxygen, and the single photo-absorption by 
0 2 and 0^. The effect of these three processes will be calculated below. 


CA37) 


The intensity at any height h is given by 

T T "A 

Iu = I e 
h o 

where I is the initial intensity and A is the optical depth. For single 


photo -absorption A is given by 


a = / ^ 

j O <•' 

where n^ are the densities of different constituents in the atmosphere and 

cr^ are the corresponding photo-absorption cross sections. 

We consider laser radiation of wavelength 2256A. The single photo- 

30 

absorption cross sections for 0 2 and C> 3 are given by 

04 *? 1 Q O 

a(0 2 ) = 4x10 cm , a(0 3 ) = 2.75x10“ cm 
Using the density profile for 0 2 and 0^ (Ref. 29) and the above cross sections 
the ratio I^/I 0 can be calculated. This ratio is plotted versus the altitude 
in Figure 6. It is seen that below 50 Km the laser beam in substantially 
absorbed by the ozone. 

In calculating the two photon absorption by 0, it should be kept in 
mind that part of the radiation absorbed are given back to the field through 
the two-photon stimulated emission. Stimulated emission should also be 
taken into account for single photo-absorption, but due to the weakness of 
the intensity the effect can be neglected. It can then be shown that 


A = 1 f c£\ T ft (* 39 ) 

J o ^ 05 " J o 

where P = N 2 /N q is a small number (Cf. Fig. 3), and has been neglected, n 
is the density of atomic oxygen. 

29 

Using values for column density of atomic oxygen , the ratio Ij /I q 
for absorption of atomic oxygen is plotted in Figure 6. It is seen that 
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A 


in the range 120-20 Km this absorption is less than 2 percent and is 
negligible compared to absorption by 0^ and 0 2> 

In Table IV distances in horizontal directions equivalent to unit 
optical depth for various altitudes due to absorption of A=2256A radiation 
by ozone are given. 
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TABLE IV. Temperature, column densities for 0, 0 and 0_, the total column 
density, the collision time T , and distances equivalent to unit optical depth 
for absorption of A=2256A radiation by ozone, X Q , as functions of the altitude. 

N(cm“ 2 ) 


h(Km) 

T(°K) 

0 

°2 

°3 

Total 

T c (sec) 

X o (Km) 

120 

360 

1.68+17 

3.92+16 

1.18+8 

6.01+17 



115 

300 

2.25+17 

7.26+16 

4.93+8 

1.23+18 

1.06-3 

4.85+9 

110 

240 

3.21+17 

1.56+17 

2.06+9 

2.50+18 

5.79-4 

1.16+9 

105 

209 

4.64+17 

3.90+17 

8.58+9 

5.11+18 

3.03-4 

2.79+8 

100 

195 

6.-45+17 

1.06+18 

3.58+10 

1.04+19 

1.55-4 

6.69+7 

95 

188 

8.73+17 

2.91+38 

1.49+11 

2.12+19 

7.76-5 

1.60+7 

90 

187 

1.04+18 

7.72+18 

6.23+11 

4.32+19 

3.79-5 

3 . 84+6 

85 

189 

1.12+18 

1.96+19 

2.60+12 

9.35+19 

1.66-5 

9.20+5 

80 

199 

1.15+18 

4.72+19 

1.09+13 

2.25+20 

6.17-6 

2.20+5 

75 

208 

1.16+18 

1.08+20 

4.52+13 

5.14+20 

2.74-6 

5.29+4 

70 

220 

1.17+18 

2.36+20 

1.89+14 

1.13+21 

1.26-6 

1.27+4 

65 

233 

1.17+18 

4.94+20 

7.96+14 

2.36+21 

6.09-7 

2.99+3 

60 

247 

1.18+18 

9.92+20 

2.97+15 

4.73+21 

3.07-7 

8.37+2 

55 

261 

1.18+18 

1.92+21 

9.29+15 

9.16+21 

1.60-7 

2.88+2 

50 

271 

1.18+18 

3.61+21 

2.86+16 

1.72+22 

8.64-8 

9.41+1 

45 

264 

1.19+18 

6.76+21 

8.94+16 

3.22+22 

4.69-8 

2.99+1 

40 

250 

1.19+18 

1.30+22 

2.74+17 

6.20+22 

2.44-8 

9.87+1 

35 

237 

1.19 + 18 

2.60+22 

7.41+17 

1.24+23 

1.20-8 

3.88 

30 

227 

1.19+15 

5.38+22 

1.68+18 

2.57+23 

5.72-9 

1.94 

25 

222 

1.19+18 

1.15+23 

3.34+18 

5.46+23 

2.66-9 

1.10 

20 

217 

1.19+18 

2.48+23 

5.61+18 

1.18+23 

1.22-9 

8.00-1 


- 39 - 


***** 4 . 


FIGURE CAPTIONS 


Figure 1. Fraction of atomic oxygen in the excited 2p 3 (^S°)3p 3 P and ionized 

states, N„/N and N,/N , as functions of the laser pulse duration x using 
2 o 3 o 

C2*i 2 

the relationship R 2r - J xl^ ' ^he l aser pulse has the characteristic 
2 

e/A=9.82 joule/m , where e is the laser's energy, and A is its cross sectional 
area. The calculated values are valid for x<x p shown by the solid curves. 

Tp is the inverse of the Doppler broadening Vp and is taken to be 1.27x10 ^sec 
corresponding to 300°K. The position of Tp is shown by the vertical solid line 

( 7 ) 2 

Figure 2. Similar to Figure 1, but for R 2r = Xpl^ * calculateci 

values are valid for x>Xp. 

Figure 3. Combination of Figures 1 and 2 for N 2 /N q values. N 2 /N q '*' s also 
the fraction of atoms at the end of each pulse that undergo spontaneous 
emission. 

Figure 4. Combination of Figures 1 and 2 for N 3 /N q values. N 3 /N 0 is als0 
the fraction of atoms that are ionized. 

3 

Figure 5. Collision time of the excited state (3p P) of oxygen with the 
atmospheric gas as functions of altitude using gas kinetic cross sections. 

Figure 6. Atmospheric absorption for A=2256A radiation due to 0^, 0 2> and 0 
as functions of the altitude for rays that propagate vertically downward. I 
and 1^ are intensities at 120 and at the height h, respectively. Absorption 
by 0 is due to two-photon processes. 
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